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I.  Introduction 

High  energy  laser  beams  propagating  through  the  atmosphere  can  be  severely 
defocused  or  deflected  by  thermal  blooming  [l].The  thermal  blooming  process  is  driven 
by  a  small  fraction  of  the  laser  energy  that  is  absorbed  by  the  molecular  and  aerosol 
constituents  of  air  [2,3]-  The  absorbed  energy  locally  heats  the  air  and  leads  to  a  decrease 
in  the  air  density  which  modifies  the  refractive  index.  The  refractive  index  variation 
leads  to  a  defocusing  or  deflection  of  the  laser  beam.  In  the  presence  of  a  transverse  wind, 
the  region  of  heated  air  is  convected  out  of  the  beam  path  and  a  steady-state  situation  is 
realized  [1].  hi  general,  however,  the  intensity  of  a  beam  undergoing  thermal  blooming 
is  a  function  of  both  time  and  spatial  position,  particularly  in  a  stagnation  zone,  where  the 
effective  wind  velocity  is  zero. 

Stagnation  zones  are  particularly  detrimental  to  HEL  propagation  since,  without 
an  effective  clearing  mechanism  for  the  heated  air,  the  strength  of  the  thermal  lens  grows 
in  time.  In  this  situation,  the  defocusing  process  is  eventually  limited  by  thermal 
conduction  or  buoyancy.  However,  by  the  time  these  processes  become  effective,  the 
beam  may  have  already  been  severely  degraded. 

Thermal  blooming  in  the  presence  of  a  stagnation  zone  was  experimentally 
observed  in  a  laboratory  experiment  using  a  10W,  CCb  laser  passed  through  an  absorption 
cell  containing  CO2  gas.  The  cell  was  pivoted  to  simulate  a  wind  profile  containing  a 
stagnation  zone.  This  experiment  was  also  modeled  using  a  code  that  solved  time- 
dependent  thermal  blooming  equations  in  the  isobaric  regime  [4], 

In  this  study  we  examine  the  effect  of  a  stagnation  zone  on  a  HEL  beam  focused 
onto  a  remote  target  in  a  maritime  atmosphere.  We  first  analyze  the  effect  of  a  stagnation 
zone  on  the  propagation  of  a  focused  beam  and  calculate  an  approximate  expression  for 
the  relative  intensity  (Strehl  ratio)  on  target  as  a  function  of  focusing  geometry  and 
stagnation  zone  position.  We  then  utilize  the  HELCAP  code  [3],  a  3D,  time-dependent 
numerical  simulation,  to  examine  the  propagation  of  a  megawatt-class  HEL  beam 
through  a  maritime  atmosphere  which  contains  a  stationary  stagnation  zone  along  the 
laser  path.  The  laser  target  configuration  used  in  these  simulations  is  shown  in  Fig.  1. 

Our  simulation  model  includes  time-dependent  thermal  blooming  as  well  as  aerosol  and 
molecular  absorption  [2,3],  turbulence  [5],  aerosol  and  molecular  scattering  [6],  and 
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aerosol  heating  and  vaporization  [7-10].  We  use  the  simulations  to  calculate  the  laser 
power  delivered  to  a  remote  target  for  different  wind  profiles  which  place  the  stagnation 
zone  at  different  locations  along  the  laser  propagation  path. 


II.  Analysis  of  HEL  Propagation  Through  a  Stagnation  Zone 

The  effect  of  a  stagnation  zone  on  HEL  propagation  in  the  atmosphere  is  highly 
dependent  on  laser  and  atmospheric  parameters  as  well  as  on  the  propagation  geometry. 
In  general,  the  problem  requires  a  numerical  simulation  to  determine  the  amount  of 
power  reaching  a  remote  target.  Before  presenting  results  of  full-scale  propagation 
simulations,  however,  it  is  useful  to  analyze  the  following  simpler  problem. 

Consider  a  focused  laser  beam  of  wavelength,  A ,  propagating  from  z  =  0  to  a 
remote  target  at  a  range  z  =  L ,  where  z  denotes  the  propagation  coordinate.  The  focal 


length  of  the  laser  beam,  Lf ,  is  such  that  Lf  >  L .  The  wind  speed,  V(z)  is  only  a 
function  of  z  and  has  a  stagnation  point  defined  by  V(z0 )  =  0 ,  where  0  <  z0  <  L.  The 


atmosphere  is  characterized  by  an  absorption  coefficient,  a  ,  and  a  scattering  coefficient. 
/? .  Turbulence  is  neglected  in  this  analysis.  This  configuration  allows  an  approximate 


expression  for  the  Strehl  ratio  (relative  intensity)  at  the  target  to  be  calculated  analytically 
as  a  function  of  both  the  location  of  the  absorbing  region  and  the  laser  focal  length. 

The  propagation  of  the  laser  beam  is  described  by  the  paraxial  wave  equation 


8  A 
~dz 


r 


i 


Vi  +  i  ^8nw  -  l 

\2  K  c  2 


A, 


(1) 


where  the  laser  electric  field  is  E  =  A  exp[  i(co0  z!  c  -  co0  /)]  e^.  / 2  +  c.c.,  co0  -2nd  A , 


ex  is  a  unit  polarization  vector  in  the  x  direction,  k0  =  n0  co0  /  c ,  n0  is  the  ambient 


refractive  index,  Vx  is  the  gradient  operator  in  the  transverse  direction,  c  is  the  speed  of 
light  in  vacuum,  Snw  is  the  change  in  refractive  index  induced  by  the  laser,  and  c.c. 
denotes  the  complex  conjugate.  In  the  absence  of  a  transverse  wind  or  slew,  there  is  axial 
symmetry  about  the  z-axis  and  all  spatially  varying  quantities  are  functions  of  ( r,z,t ) 
only,  where  r  is  the  radial  coordinate  and  t  is  time.  Writing 

A(r,z,t )  =  2?(r,z,f)exp[/0(r,z,/)],  where  B ,  6,  and  8 nw  are  assumed  to  be  real,  leads 
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to  the  following  equations  for  the  laser  intensity  1  (~  B1)  and  the  transverse 


wavenumber  kL  =V±0: 


dl_ 

dz 


■vAki)-pi, 


(2) 


dkL 

dz 


8  K 


vi-i-h/vi/  -  (v,/)2)! 


-7“ki-Vik1  +  —VLSnn. 
kn  c 


(3) 


On  the  right  side  of  Eq.  (3),  the  first  term  is  of  order  ~  A  /  R}  where  R  is  the  laser  spot 


size.  The  second  term  is  of  order  ~  R/(L2fA)  for  a  focused  beam  and  the  third  term  is  of 


order  ~  5  nTB  /{A  R ) .  We  consider  parameters  such  that  R  ~  50  cm ,  A  ~  1  pm ,  and 
Lf  ~  5  km  .  In  this  parameter  regime,  the  second  term  is  much  larger  than  the  first  term. 

The  third  term  is  much  larger  than  the  second  term  for  5  nw  >  10“8 ,  which  is  typically 
encountered  within  a  stagnation  zone.  Hence,  for  the  parameter  regime  of  interest,  we 
neglect  the  first  and  second  terms  on  the  right  side  of  Eq.  (3)  and  integrate  to  solve  for 
k_L ,  i-e.. 


kx(r,z,0  *  kx0(r)  +  ^  \v xSnTB{r,z' ,t)dz' ,  (4) 

C  0 

where  ki0(/')  is  the  initial  transverse  wavenumber  at  z  =  0.  Equation  (4)  represents  kx 

in  the  geometric  optics  limit. 

Equation  (2)  can  be  integrated  to  yield 

I(r,z,t)  =  exp]  -  —  J—  1  -Vx-[kxe r,z\t)I{r,z\t)]dz\ ,  (5) 

[  w0  0J/(r,z',O  J 

which  is  an  integral  equation  for  I(r,z,t)  that  can  be  solved  approximately  by  replacing 

I(r,z,t )  in  the  right  side  by  I0(r,z),  i.e.,  the  unperturbed  intensity  for  a  homogeneous 

medium  (S  nw  =0).  The  unperturbed  intensity  satisfies 

I0(r,z)  =  7(r,0)e  ^exp J-—  J— V,  \k^r)I0{r,z')]dz'  .  (6) 

Substituting  Eq.  (4)  into  Eq.  (5)  yields  the  relative  intensity  [11] 
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(7) 


The  refractive  index  change  due  to  thermal  blooming  is  given  by 
5n TB  =  (na  - 1  )Sp  /  p0 ,  where  p0  is  the  mass  density,  and  Sp  is  the  perturbed  mass 
density.  In  the  isobaric  regime,  the  mass  density  change  is  given  by  [1,3] 


CnT  —  +  V-V  - 
p  0  a* 


- - V2  5p  =  -aw(r,t)I , 

CpPo  J 


(8) 


where  k  is  the  thermal  conductivity,  V  is  the  effective  wind  velocity,  Cp  is  the  specific 
heat  at  constant  pressure,  T0  is  the  ambient  temperature,  and  1  is  the  time-averaged  laser 
intensity.  In  general,  the  absorption  coefficient,  aTB ,  contains  both  aerosol  and  molecular 
contributions  and  can  be  time  dependent.  For  the  simplified  analysis  of  this  Section, 
however,  we  take  it  to  be  constant  in  time.  The  isobaric  approximation  is  valid  for  times 
greater  than  the  hydrodynamic  time  R / Cs,  where  C\  is  the  acoustic  speed.  We  consider 

time  scales  less  than  the  thermal  conduction  time,  Cpp0R 2  Ik,  and  neglect  the  thermal 
conduction  term  in  Eq.  (8). 

For  the  purpose  of  obtaining  an  approximate  closed  form  expression  for  Eq.  (7), 
we  solve  Eq.  (8)  approximately  by  writing  the  operator  V-V  »  1/  xc (z)  where 

tc(z )  =  R(z)/\  V(z)  |  is  the  clearing  time.  Using  this  approximation,  we  can  integrate  Eq. 
(8)  to  obtain 


Snw(r,z,t)  *  -y7BaI0(r,z)rc(z){\  -  exp[- 1 / rc (z)]} , 


(9) 


where  yTB  =  (n0  -\)/(CpT0p0)  . 


The  unperturbed  intensity  is  taken  to  have  a  Gaussian  transverse  profile 


(10) 
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where  700  denotes  the  intensity  on  axis  at  z  =  0  .  The  spot  size  R  is  a  function  of  z  and 
the  initial  spot  size  is  denoted  by  R0  =  R(z  =  0) .  Substituting  Eqs.  (9)  and  (10)  into  Eq. 
(7),  results  in  an  expression  for  the  on-axis  Strehl  ratio  at  z  =  L ,  i.e.. 


I(r  =  0,L,t) 
h(r  =  0,L) 


exp 


L, 

- %IwRlyna\dz'\dz"GV',t ) 


0  0 


(11) 


where  G(z,t)  =  rc(z)[l  -  exp(-//'rt(z))]exp(-/?z)/./?4(z).  To  obtain  an  approximate 


analytic  expression  for  the  right  side  of  Eq.  (1 1),  we  note  that  for  a  weakly  focused  beam, 
the  integrand  G  is  sharply  peaked  and  grows  linearly  with  time  in  the  vicinity  of  the 
stagnation  zone,  i.e.,  G(z0,t )  =  texp(-/? z0)/  R4(z0) .  Hence,  we  can  approximate  the 


integration  over  z"  in  Eq.  (1 1)  by 

\dz"  G(z",t )  »  G(zo,OAz0(z’-zo),  (12) 

o 

where  Az  is  the  characteristic  width  of  the  stagnation  region  and  0  is  the  Heaviside  step 
function. 

Using  Eqs.  (12)  in  Eq.  (1 1)  results  in 


I(r  =  0,L,O 
h(r  =  0,L) 


8ywaI00L2t  { Az^  _  ^  R\  exp(-/?z0) 

’  L 


R 


0 


L 


y\ 


j 


R4(z  o) 


(13) 


The  width  of  the  stagnation  region,  Az ,  is  in  general  a  function  of  time  and  also 
the  spatial  variation  of  the  wind  profile.  We  can  define  the  stagnation  zone  as  the  region 
where  the  condition  t  /  rc  (z)  «  1  is  satisfied.  In  this  case,  the  width  of  the  stagnation 


zone  can  an  be  defined  according  to  t  /  rc(z0  +  Az)  =  1 . 

We  now  consider  specific  analytic  expressions  for  the  wind  profile  V(z)  and  laser 
spot  size  R(z) .  For  a  slewed  laser  beam,  the  wind  profile  can  be  written  as 

V(z)  =  Vmni/(l~z/z0),  (14) 


where  Vwmd  is  the  ambient  wind  velocity,  and  the  laser  slew  rate  ,  6  is  related  to  z0 
according  to  0  =  Vmnd  /  z0 .  Note  that  for  this  model,  a  stagnation  zone  at  the  transmitter 
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(  z  =  0 )  implies  an  infinite  slew  rate.  Assuming  that  the  laser  spot  size  does  not  vary 
appreciably  within  the  stagnation  zone. 


Az  *  z0R(za)/(Vwmdt ) 


(15) 


for  the  wind  profile  of  Eq.  (14). 

For  the  linearly  varying  wind  profile  described  above,  Eq.  (13)  can  be  written  as 


(16) 


where 


(17) 


describes  how  the  location  of  the  stagnation  zone  affects  the  relative  intensity  on  target. 
The  larger  the  value  of  /(z0) ,  the  lower  the  intensity  on  the  target. 

We  can  derive  an  approximate  functional  form  for  R(z )  to  be  used  in  Eq.  (17). 
Substituting  Eq.  (10)  into  Eq.  (6),  the  laser  spot  size  can  be  shown  to  satisfy 
R(z)  =  /?0[1  -  z(l  -  r/)/  L]  for  z  «  Lf  ,  where  the  parameter  77  s  R(L)I  R0,  i.e.,  the 

ratio  of  the  spot  size  on  target  to  the  initial  spot  size.  In  terms  of  the  focal  length, 

77  =  1  -  (LI  Lf) .  Phenomenologically,  77  can  also  parameterize  focusing  limitations  due 

to  beam  quality  or  atmospheric  turbulence. 

Figure  2  plots  the  distortion  function  /  versus  z0  for  various  values  of  77 .  We 
assume  a  normalized  scattering  coefficient  of  (3  L  =  0.25  .  For  a  beam  that  is  focused 
near  the  target  plane,  the  distortion  at  the  target  plane  is  largest  when  the  stagnation  zone 
is  close  to  the  target.  In  contrast,  when  the  beam  is  collimated  (77  ~  1 ),  the  distortion  at 
the  target  plane  is  greater  when  the  stagnation  zone  is  closer  to  the  laser  source.  Note  that 
the  case  z0  =  0 ,  which  corresponds  to  an  infinite  slew  rate  given  the  assumed  wind 

profile  of  our  model,  is  not  physically  realizable. 

III.  Numerical  Simulation 

The  propagation  code  used  in  this  study  is  HELCAP  (High  Energy  Laser  Code  for 
Atmospheric  Propagation),  which  is  a  fully  time-dependent,  3D  code  developed  at  the 
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Naval  Research  Laboratory.  HELCAP  models  the  propagation  of  continuous  and  pulsed 
HEL  beams  through  various  atmospheric  environments.  It  includes  the  effects  of  aerosol 
and  molecular  scattering,  aerosol  heating  and  vaporization,  thermal  blooming  due  to  both 
aerosol  and  molecular  absorption,  and  atmospheric  turbulence.  A  more  detailed 
description  of  the  code  can  be  found  in  Refs.  [3]  and  [12]. 

HELCAP  solves  a  nonlinear  Schrodinger-like  equation  for  the  complex  laser 
electric  envelope,  A(x,y,z,t),  which  is  of  the  form 


d_A 
d  z 


ic. 

~2co 


V{A 


\  i  —  {Sn  , 
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7  +  8nIB)  +  P) 


1^7 


(12) 


where  A(x,y,z,t )  is  the  complex  laser  electric  field  amplitude  and  the  terms  denoted  by 
^Sj  represent  other  physical  processes  that  are  included  in  the  code  but  are  not 

j 

important  for  the  parameter  regime  considered  here.  The  laser  intensity  is 
I  =  c  A  A*  / 8 x .  The  absorption  and  scattering  coefficients  aTB  and  /? ,  respectively, 
contain  both  molecular  and  aerosol  contributions.  In  general,  they  are  time  dependent  due 
to  aerosol  heating  and  vaporization  by  the  laser  as  described  in  Ref.  [3].  The  quantities 
dnT  and  SnrB  denote  the  refractive  index  variation  due  to  atmospheric  turbulence  and 
thermal  blooming  respectively.  Atmospheric  turbulence  is  modeled  in  the  usual  manner 
using  phase  screens  for  which  the  scale  sizes  of  the  index  fluctuations  are  described  by  a 
Kolmogorov  spectrum  characterized  by  the  parameter  Cn2  [5],  The  mass  density  change 
due  to  thermal  blooming  is  solved  for  according  to  Eq.  (8),  where  the  total  absorption 
coefficient ,  aTB ,  is  time-dependent  and  contains  both  aerosol  and  molecular 
contributions  in  the  manner  described  in  Ref.  [3]. 

For  the  numerical  simulations  presented  here,  we  consider  a  slewed  HEL  beam 
focused  at  a  range  of  L  =  5  km  in  the  presence  of  a  uniform  transverse  wind.  The  z-axis 
defines  the  direction  of  propagation  and  the  wind  and  slew  velocities  are  taken  to  be 
along  the  y-axis.  The  simulation  geometry  is  shown  in  Fig.  1.  The  high  energy  laser  beam 
has  an  initial  field  profile  given  by  A  =  AQf  (r)g-(/)exp(-r2  /  R* ) ,  where 

f(r)  =  exp[-(2r/  D)f  ]  ,£  =  20,  limits  the  transverse  extent  of  the  beam  to  the  aperture 
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diameter,  D.  The  function  g(t )  is  the  initial  temporal  profile  of  the  beam  which  contains 
a  pointing  jitter  characterized  by  an  angular  spread  of  2  grad  and  a  white  noise  temporal 
spectrum.  The  transmitted  power  at  the  source  is  denoted  by  Pr  . 

The  laser  is  propagated  through  a  maritime  environment  in  which  the  aerosol 
distribution  is  modeled  using  the  Advanced  Navy  Aerosol  Model  (ANAM)  [13]  and  the 
molecular  absorption  coefficient  is  calculated  using  MODTRAN4  [14].  The  details  of  the 
atmospheric  model  and  the  parameters  used  to  characterize  a  maritime  environment  are 
discussed  in  Ref.  [3].  The  ambient  absorption  and  scattering  coefficients  for  propagation 
at  X  =  2.141  /jm  are  an  =  6xl0'3  km'1  and  /?  =  0.05  km'1 ,  respectively.  In  the 
simulations  these  quantities  vary  in  space  and  time  due  to  aerosol  vaporization. 

The  effective  wind  speed  along  the  propagation  path  is  V(z )  =  Vmnd  -  9  z ,  where 

z  is  the  distance  from  the  laser  source,  VwtM  is  the  ambient  wind  speed,  and  9  is  the  laser 

beam  slew  rate.  For  positive  9,  i.e.,  slewing  in  the  direction  of  the  wind,  a  stagnation 
zone  is  created  at  location  z0  =  Vwind  19.  In  what  follows,  we  investigate  the  effect  of 

this  stagnation  zone  on  the  propagation  of  a  HEL  beam  for  various  slew  rates  which 
change  the  distance  from  the  transmitter  to  the  stagnation  zone. 

IV.  Results 

We  consider  a  laser  of  wavelength  X  =  2.141pm  and  power  PT  =1.5  MW 
propagating  through  a  maritime  atmosphere  with  Vwnd  =  5  m  /  sec .  The  power  delivered 

to  a  target  of  circular  area  100  cm2  at  a  range  of  5  km  is  used  as  a  figure  of  merit.  In  the 
simulations  that  follow,  the  target  is  always  located  such  that  the  peak  laser  fluence  at  5 
km  range  is  at  the  center  of  the  target.  The  average  power  reaching  the  target  is  defined 


by 
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where  the  dwell  time  vdMl  =  1  sec,  and  dxdy  is  the  differential  cross  section  which  is 
integrated  over  the  target  area.  The  total  laser  energy  reaching  the  target  is 


Figure  3  plots  ^Ptarget  j  versus  slew  rate,  9 ,  for  the  propagation  configuration  of 

Fig.  1.  A  negative  slew  denotes  that  the  slew  direction  is  opposite  to  the  wind  direction. 

In  this  case,  there  is  no  stagnation  zone  between  the  laser  and  the  target.  A  positive  slew 
denotes  slew  in  the  direction  of  the  wind.  In  this  case,  it  is  possible  to  have  a  stagnation 
zone  somewhere  along  the  propagation  path.  For  example,  a  slew  rate  of  1  mrad/sec 
places  the  stagnation  zone  at  the  target,  while  a  slew  rate  of  5  mrad/sec  places  the 
stagnation  zone  1  km  from  the  laser  source.  Figure  3  shows  that  in  the  absence  of  slew 
(i.e.,  uniform  wind),  the  power  on  target  is  )  »  0.7  MW  and  the  propagation 

efficiency  is  ~  50%.  ^P(argot  ^  increases  relative  to  tins  value  when  the  slew  is  opposite  to 

the  wind  direction  because  the  effectively  larger  wind  speed  along  the  entire  propagation 
path  mitigates  thermal  blooming.  For  a  slew  rate  9  =  -5  mrad/sec ,  for  example, 

(P{as ,,,  j  «  1.1  MW ,  which  corresponds  to  a  propagation  efficiency  of  >  70%. 

When  the  slew  is  increased  in  the  direction  of  the  wind,  the  power  on  target 
decreases  sharply  for  9  <  1  mrad/sec  and  reaches  a  minimum  when  9  =  1  mrad/sec  , 
i.e.,  when  the  stagnation  zone  is  at  the  target.  At  this  minimum,  /p|areet  ^  ~  0.07  MW , 

which  corresponds  to  a  propagation  efficiency  of  <  5%.  Increasing  the  slew  rate  further 
brings  the  stagnation  zone  closer  to  the  laser  source,  where  the  intensity  is  lower,  and 
increases  the  effective  wind  velocity  near  the  target.  The  result  is  that  the  power  on  target 

increases  with  slew  rate.  For  a  slew  rate  of  9  =  5  mrad/  sec ,  (Ptareet  ^  ~  0.7  MW  .  These 

results  are  consistent  with  the  calculation  in  Sec.  II,  i.e.,  for  a  beam  focused  onto  the 
target  plane,  a  stagnation  zone  close  to  the  target  will  result  in  the  largest  beam  distortion. 

Figure  4  plots  the  laser  fluence  on  target  calculated  over  a  1  sec  dwell  time.  Panel 
labels,  (a),  (b),  (c),  and  (d)  correspond  to  the  different  slew  rates  indicated  in  Fig.  3. 

Panel  (a)  corresponds  to  a  slew  rate  9  =  -5  mrad/  sec ,  which  gives  the  maximum  power 
on  target  over  the  range  of  slew  rates  considered.  The  beam  is  well-focused  and  exhibits 


10 


practically  no  thermal  blooming  effects.  Panel  (b)  shows  the  beam  profile  when  the  slew 
rate  is  zero,  i.e.,  a  uniform  wind  profile.  In  this  case  the  beam  fluence  profile  shows  the 
characteristic  crescent  shape  associated  with  thermal  blooming  with  the  centroid  of  the 
beam  deflected  opposite  to  the  direction  of  the  wind.  Panel  (c)  shows  the  beam  profile  on 
target  for  the  case  when  the  stagnation  zone  is  at  the  target.  The  beam  is  severely 
defocused  as  a  result  of  thermal  blooming.  The  peak  in  fluence  near  (x  =  0,  y  =  0)  is  due 
to  the  leading  edge  of  the  beam  (t  <  20  msec),  which  is  not  affected  by  the  time 
dependent  thermal  blooming  process.  Panel  (d)  shows  the  beam  profile  when  the 
stagnation  zone  is  1  km  from  the  target.  In  this  case,  the  effective  wind  velocity  near  the 
target  is  in  the  opposite  direction  relative  to  that  of  panel  (b).  The  value  of  ^Ptarget  ^  is 

similar  to  that  of  panel  (b)  but  with  the  beam  deflected  in  the  opposite  direction. 

Figure  5  shows  the  time  dependence  of  the  beam  intensity  at  the  target.  Figure  4a 
corresponds  to  a  slew  rate  of  9  =  -5  mrad/  sec ,  when  the  power  on  target  is  maximum.  In 
this  case  the  beam  is  well-focused  over  the  entire  dwell  time.  The  transverse  profile  is 
mainly  affected  by  turbulence  and  pointing  jitter.  Figure  4b  corresponds  to  a  slew  rate  of 
6  =  1  mrad / sec  ,  which  results  in  the  minimum  power  on  target.  For  this  case,  it  is  seen 
that  severe  defocusing  occurs  for  times  >  20  msec,  while  the  first  ~10  msec  of  the  beam 
is  relatively  unaffected  by  thermal  blooming. 

V.  Conclusions 

Stagnation  zones  can  significantly  degrade  the  propagation  efficiency  of  a  high 
energy  laser  (HEL)  through  the  atmosphere.  We  have  analyzed  the  propagation  of  a 
focused  HEL  beam  through  a  model  atmosphere  containing  a  stagnation  zone.  Our 
analysis  shows  that  the  amount  of  laser  degradation  caused  by  a  stagnation  zone  is  highly 
dependent  on  the  focusing  geometry  of  the  laser  beam  and  the  location  of  the  stagnation 
zone  along  the  propagation  path.  For  a  collimated  laser  beam,  a  stagnation  zone  near  the 
transmitter  causes  the  largest  distortion  of  the  laser  beam  on  the  target.  For  a  laser  that  is 
focused  onto  the  target,  a  stagnation  zone  near  the  target  causes  the  largest  laser 
distortion  and  the  fastest  decrease  in  laser  intensity.  This  finding  is  in  agreement  with  the 
experimental  results  of  Ref.  [4]. 
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We  use  a  full  scale  numerical  simulation  (HELCAP)  to  model  the  propagation  of 
a  megawatt-class  HEL  beam  through  a  realistic  maritime  atmosphere  containing  a 
stationary  stagnation  zone.  The  HEL  beam  is  focused  onto  a  remote  target.  The 
stagnation  zone  is  created  by  slewing  the  HEL  beam  in  the  direction  of  the  wind  and 
location  of  the  stagnation  zone  is  varied  by  changing  the  slew  rate.  Consistent  with  our 
analysis,  the  simulations  show  that  the  power  on  target  is  minimized  when  the  slew  rate  is 
such  that  the  stagnation  zone  is  located  near  the  target  plane.  In  this  case,  propagation 
efficiency  can  be  reduced  by  an  order  of  magnitude  relative  to  an  unslewed  beam. 
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Figure  1 :  Schematic  diagram  showing  the  simulation  geometry.  A  HEL  beam  of  power 
Pt  is  focused  at  a  range  of  L  =  5  km  in  the  presence  of  a  uniform  wind  with  speed 
Vw  =  5  m/sec  .  The  beam  is  slewed  in  the  direction  of  the  wind  to  create  a  stagnation 
zone  along  the  propagation  path. 
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Figure  2:  Distortion  function  /  as  defined  by  Eq.  ()  plotted  versus  distance  to  stagnation 
zone  z0  for  various  values  of  the  focusing  parameter  rj  -  R(L)/ R0  where  R(L)  is  the 
unperturbed  laser  spot  size  at  the  target  plane  ( z  =  L ),  R0  is  the  initial  spot  size  at  z  =  0  , 
and  /3L  =  0.25 . 
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Figure  3:  Average  power  on  target,  ,  versus  slew  rate  for  k  =  2.141pm, 

Pr  =  1.5  MW  ,  Vwlnd  =  5  m/sec,  tM,  =  1  sec,  and  L  =  5  km  .  Points  (a),  (b),  (c),  and 
(d)  are  labeled  to  correspond  with  fluence  distributions  shown  in  Fig.  3. 
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Figure  4:  Laser  fluence  profiles  at  the  target  range  for  slew  rates  (a)  6  =  -5  mrad/  sec  , 
(b)  0  =  0 ,  (c)  6  =  1  mrad/  sec,  and  (d)  0  =  5  mrad /sec  .  Parameters  are  the  same  as 
those  of  Fig.  2.  Fluence  is  calculated  over  a  1  sec  dwell  time. 
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Figure  5:  Laser  intensity  in  the  y=0  plane  versus  transverse  coordinate,  x,  and  time,  t,  for 
slew  rates  (a)  6  =  -  5  mrad/sec  and  (b)  6  =  1  mrad/sec .  Parameters  correspond  to 
those  of  Fig.  2. 
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